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Abstract. The volume and surface effects in the nuclear local energy density and the volume and sur- 
face components of the pairing interaction are discussed in the context of the mean-field, Hartree-Fock- 
Bogoliubov description of atomic nuclei. Predictions of properties of exotic nuclei close to the particle drip 
lines are presented. 

PACS. 21.60.Jz Hartree-Fock and random-phase approximations - 21.10.Dr Binding energies and masses 
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1 Introduction 

The mean-field methods are very successful in describing 
and predicting properties of nuclei across the chart of the 
nuclides. This is especially true for heavy nuclei, where 
the bulk properties of nuclear matter dominate over the 
surface effects. However, when details of nuclear structure 
are considered, a correct description of the nuclear surface 
is essential. Moreover, the surface region may give us in- 
valuable information on the nature and strength of nuclear 
effective interactions in channels that are inaccessible by 
considering infinite systems (i.e., nuclear matter). In the 
present report, we briefly discuss several aspects of the 
surface effects in the particle-hole (Sect. ||) and particle- 
particle channels (Sect. ||), as well as the deformation ef- 
fects (Sect. |). 



2 Volume and surface components of the 
energy density 

Without any detailed microscopic knowledge of the nu- 
clear effective interactions, we can rely on general prop- 
erties of saturating fermion systems to assume that the 
total energy £ of a nucleus is an integral of a local energy 
density Tt(r). Such a conjecture is a basis for the so-called 
Local Density Approximation (LDA) that has been exten- 
sively used in the context of atomic and molecular physics. 
In nuclear physics LDA has been employed in the form of 
the Skyrme-HF approximation, in which the total energy 
£ is given as 



or equivalently 
£= I d 3 r 



Wkin(f) + Hskyrmc(r) + Ws-o(^) 
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The kinetic (7Y k in), Skyrme (Wskyrme), spin-orbit (Hs-o), 
Coulomb (Hcoui), and pairing (7i pa ir) densities are func- 
tions of several local densities: 
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where, e.g., po=p n +p p and pi=p n —p P are the isoscalar 
and isovector particle densities, respectively, and p n and 
p p are the corresponding neutron and proton densities. 
For complete definitions of other densities and coupling 
constants appearing in expressions (0)~(@), the reader is 
referred, e.g., to Refs. |,|. (In Eq. (g) we have omitted 
the term depending on the tensor spin-current density be- 
cause below we use only the Skyrme force in which this 
particular term was neglected.) 

In the lower panel of Fig. [y the above five energy den- 
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(1) sities are plotted for 12U Sn, together with their sum Ht 
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Fig. 1. Total energy density Htot (full circles) calculated 
within the Skyrme-HF-SLy4 method for 120 Sn, and plotted 
together with its five components given by decompositions de- 
fined in Eqs. (g) (lower panel) and ffl (upper panel). 



Calculations were performed for the Skyrme interaction 
SLy4 H and the volume pairing force. (See Ref. ||) for 
details of the calculations.) The presented results are very 
generic, and identical qualitative results are obtained for 
any other nucleus or interaction. 

Only three out of five components significantly con- 
tribute to the total energy density, namely, the kinetic, 
Skyrme, and Coulomb densities. The remaining two have, 
of course, a decisive influence on detailed properties of 
nuclei, however, they are almost invisible in the scale of 
Fig. and for the sake of the following discussion can be 
safely put aside. We can also see that the kinetic energy 
density dominates in the surface region; indeed, both the 
Skyrme and Coulomb terms simultaneously go to zero at 
a distance that is by about 1 fm smaller than the place 
where the kinetic energy vanishes. 

Therefore, in order to analyze the energy relations at 
the nuclear surface, it is essential to consider surface prop- 
erties of the kinetic energy density. Semiclassical methods 
are not appropriate to separate the volume and surface 
contributions to 7ikin, because such approaches are not 
valid beyond the classical turning point. Therefore, one is 
often fitting the volume and surface terms to reproduce 
the microscopic values of Ttkin (see, e.g., ||). From such 
analyses it turns out that the volume contribution is very 
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Fig. 2. Microscopic kinetic energy density T^kin (full cir- 
cles) compared to the nuclear-matter (volume) approximation 
^kin-nm (open circles). 



well described by the nuclear-matter expression 



4/3 

r nm_^ (o ^5/3 



(8) 



Hence, in the present study we simply consider the re- 
maining part of the kinetic energy density to be the sur- 
face contribution. 

In Fig. |^ we compare the microscopic kinetic energy 
density 7ikm, Eq. (|3), with the corresponding nuclear- 
matter (volume) contribution Hkin-nm obtained by replac- 
ing r = r n + t p with r$ m = T° m + r° m . The difference 
between the two curves gives our surface contribution to 
the kinetic energy density. 

Based on these arguments, we can now rearrange terms 
in £ in such a way as to single out the volume and surface 
contributions: 



Wvoi(r) +H suli (r) +H S -o(r) 
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and H vo \(r) + H snT t(r) = H k i n (r) + Wskyrme(r). In the 
Skyrme energy density (Q) the effective-mass terms have 
been separated into the volume and surface parts accord- 
ing to the prescription defined above for the bare-mass 
terms. 
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In Fig. |l| (upper panel) are plotted the volume ( [To| ) 
and surface ( |ll| ) contributions to the 120 Sn total density 
energy. It is clear that in the surface region of this nucleus 
(5-7 fm) these two contributions are of a similar magni- 
tude and opposite sign. Therefore, the nuclear surface can- 
not simply be regarded as a layer of nuclear matter at low 
density. In this zone the gradient terms (absent in the 
nuclear matter) are as important in defining the energy 
relations as those depending on the local density. 

This observation exemplifies the difficulties in extract- 
ing the pairing properties of finite nuclei from the nuclear- 
matter calculations. In particular, the nuclear-matter pair- 
ing intensity, calculated at densities below the saturation 
point, need not be the same as the analogous intensity at 
the surface of a nucleus. The nuclear-matter and neutron- 
matter calculations of the pairing gap (see Ref. || for a 
review) performed at various densities by using very ad- 
vanced and sophisticated methods, as well the best bare 
NN forces, can therefore be, at most, considered as weak 
indications of what might be the actual situation in nu- 
clei. In this respect, calculations in semi-infinite matter 
that recently became available may provide much more 
reliable information. 
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Fig. 3. Radial strength factor / pa ; r of the density-dependent 
delta interaction, Eq. ([l4|), as a function of p=po for several 
values of a. The value of p c was assumed to be 0.16 fm~ 3 . At 
each value of a, the strength Vb was adjusted to reproduce 
the neutron pairing gap in 120 Sn. The inset shows / P air/|Vb| 
as a function of dimensionless normalized density po / p c (from 
Ref. H). 



3 Volume and surface pairing interactions 

Without having at our disposal microscopic first-principle 
effective pairing interactions (with surface effects included 
as discussed in Sec. |^), one uses in the particle-particle (p- 
p) channel a phenomenological density-dependent contact 
interaction. As discussed in a number of papers (see, e.g., 
Refs. HHII), tne presence of the density dependence in 
the pairing channel has consequences for the spatial prop- 
erties of pairing densities and fields. The commonly used 
density-independent contact delta interaction, 

V v s ol (r,r') = V S(r-r'), (12) 

leads to volume pairing. A simple modification of that 
force is the density-dependent delta interaction (DDDI) 

y/ urf (r,r') = / P ah.(r)<5(r-r'), (13) 
where the pairing-strength factor is 

f pai r(r)=V {l-[p (r)/p c ] a } (14) 

and Vb, p c , and a are constants. If p c is chosen such that 
it is close to the saturation density, p c fapo(r—0), both 
the resulting pair density and the pairing potential are 
small in the nuclear interior, and the pairing field becomes 
surface-peaked. By varying the magnitude of the density- 
dependent term, the transition from volume pairing to 
surface pairing can be probed. A similar form of DDDI, 
also containing the density gradient term, has been used 
in Refs. |||l4[|. 

Apart from rendering the pairing weak in the interior, 
the specific functional dependence on p used in Eq. ( [Tij ) is 



not motivated by any compelling theoretical arguments or 
calculations. In particular, values of power a were chosen 
ad hoc to be either equal to 1 (based on simplicity), see, 
e.g., Refs. |l^,^6|, or equal to the power 7 of the Skyrme- 
force density dependence in the p-h channel 

The dependence of results on a was studied in Ref. [[[7) 
within the Hartree-Fock-Bogoliubov (HFB) approach. We 
considered four values of a=l, 1/2, 1/3, and 1/6 that cover 
the range of values of 7 used typically for the Skyrme 
forces. For p c we took the standard value of 0. 16 fm~ 3 , and 
the strength Vq of DDDI was adjusted according to the 
prescription given in Ref. |)|, i.e., so as to obtain in each 
case the value of 1.245 MeV for the average neutron gap 
in 120 Sn. The resulting pairing-strength factors ( |l4| ) are 
shown in Fig. || as functions of density p=po for the four 
values of the exponent a. It is seen that for p>0.04fm~ 3 
the pairing-strength factor / pa j r is almost independent of 
the power a. At low densities, however, the pairing interac- 
tion becomes strongly dependent on a and very attractive 
at p— >0. The pattern shown in Fig. || indicates that pair- 
ing forces characterized by small values of a should give 
rise to pair fields peaked at, or even beyond, the nuclear 
surface (halo region) where the nucleonic density is low. 

The main conclusion of Ref. jlTj is that, due to the self- 
consistent feedback between particle and pairing densities, 
the size of the neutron halo is indeed strongly influenced 
by pairing correlations; hence, by the pairing parametriza- 
tion assumed. Consequently, experimental studies of neu- 
tron distributions in nuclei are extremely important for 
determining the density dependence of pairing interaction 
in nuclei. At the same time, the strong low-density de- 
pendence of the pairing force, simulated by taking very 
small values of a in DDDI, is unphysical. The present ex- 
perimental data are consistent with about l/2<a<l. In 
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Fig. 4. Comparison between the experimental two-neutron separation energies S2N and neutron pairing gaps An (upper left 
panels, based on masses from Ref. llq]), and the corresponding results of the spherical HFB method for the Skyrme SLy4 force 
|J and five different versions of the zero-range pairing interaction (see text). 



this context, it is interesting to note that excellent fits to 
the data were obtained in Refs. Jl^,0 by taking a=2/3. 
However, at present there is no theoretical argument why 
the density dependence should be even taken in a form of 
the power law. 

Moreover, the pairing interaction is most likely of an 
intermediate character between the volume ( |l2| ) and sur- 
face forms (|l|). (See Refs. for recent analyses.) 
In particular, the force which is a fifty-fifty mixture of 
both types, 



^ ix (r,r') = -(^ 



vol 



vLd=V [l-^j S(r-v'), 

_ . (15) 

performs quite well p0( ] in reproducing the general mass- 
dependence of the odd-even mass sta ggering parameter 



Z\( 3 ) centered at odd particle numbers [pi, 22 

Figure || illustrates the role of using different types 
of the pairing interaction to predict the two-neutron sep- 



aration energies and neutron pairing gaps, respectively, 
in very neutron-rich isotones around N=82. The experi- 
mental values were calculated based on the interim 2001 
evaluation of atomic masses . 

Figure 4 nicely illustrates the effect of the so-called 
shell quenching in heavy nuclei p3[ , i.e., the vanishing of 
the effective distance between the neutron single-particle 
levels above and below a magic neutron number when ap- 
proaching the neutron drip line. The difference between 
the two-neutron separation energies above and below N = 
82 very well visualizes this effect. In fact, the experimen- 
tal data show an apparent opposite effect; however, this 
is caused by the fact that the data are available only for 
Z>50. When approaching the magic proton number, the 
neutron magic gap is slightly enhanced pi) . This effect 
is entirely absent in calculations that do not include any 
effects of correlations and deformations. 

Nevertheless, for Z<50 the effect of the shell quench- 
ing is very well visible in the calculations. Moreover, the 
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magnitude of the effect is very strongly influenced by the 
type of pairing force used. For the volume pairing force 
(|l2|), the effect is rather weak and the magic gap V=82 
is still visible even at the very drip line. However, for the 
surface pairing force (|l^) the shell gap goes to zero much 
earlier, and this tendency is accentuated for pairing forces 
that are stronger at small densities (for smaller powers of 

For the neutron pairing gap (Fig. 0) the experimen- 
tal data that exist for Z>50 do not indicate any definite 
change in the neutron pairing intensity with varying pro- 
ton numbers. However, the surface pairing interactions 
(bottom panels) give a slow dependence for Z>50 that 
is dramatically accelerated after crossing the shell gap at 
Z=50. On the other hand, the volume and intermediate- 
type pairing forces predict a slow dependence all the way 
through to very near the neutron drip line. It is clear that 
measurements of only several nuclear masses for Z<50 
will allow us to strongly discriminate between the pairing 
interactions that have different space and density depen- 
dencies. 



4 Deformation of drip-line nuclei 

Deformation of nuclei near the drip lines is a difficult and 
open problem in nuclear structure physics. It requires a si- 
multaneous description of particle-hole, pairing, and con- 
tinuum effects — the challenge that only very recently 
can be addressed by mean-field methods. Deformability 
of nuclei plays a decisive role in determining particle sep- 
aration energies and decay rates, and hence is crucial for a 
description of nuclear processes in a stellar envi ronme nt 
Very recently we have developed methods 
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to 

approach this problem by using the local-scaling point 
transformation that allows us to modify asymptotic prop- 
erties of the deformed harmonic oscillator wave functions. 
The resulting single-particle bases are very well suited for 
solving the HFB equations for deformed drip-line nuclei. 

Calculations of a complete HFB mass table are now in 
progress and will be reported in separate publications p7| , 
P8| . Here we only show a sample of the results obtained 
for the SLy4 Skyrme interaction and the intermediate- 
type pairing force (|T^) . A rather restricted size of the har- 
monic oscillator basis, limited to AT s h=14 spherical shells, 
was used, while the continuum states were included up to 
60 MeV. Figure ^ shows the two-neutron separation ener- 
gies (top panel) , the neutron pairing gaps (middle panel) , 
and the deformations j3 (bottom panel) for 1553 particle- 
bound even-even nuclei with Z<108 and iV<188. 

At the two-neutron drip line, one can see a very inter- 
esting effect of negative two-neutron separation energies 
for particle-bound (negative Fermi energy) nuclei. This is 
the result of a sudden change in configuration when ap- 
proaching the drip line, that is caused by the fact that the 
ground-state configuration may become particle-unbound 
earlier than the excited one [g|[|2^]. A similar effect oc- 
curs also in the heavy proton drip-line nuclei, where se- 
quences of oblate ground states are obtained. Another ef- 
fect at the proton drip line is related to long sequences 



of proton-magic (e.g., Z=50 and 82) isotopes intruding in 
the territory of unbound nuclei. This is the result of the 
vanishing pairing correlations, for which the proton Fermi 
energy coincides with the last occupied level, while in the 
neighboring nuclei it is located higher. 
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